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1. Introduction 

A sophisticated analysis of EGRET data (Dixon et al. 1998) has found evidence 
for gamma-ray emission from the galactic halo. Filtering the data with a wavelet 
expansion, Dixon et al. have subtracted an isotropic extra-galactic component and 
expected contributions from cosmic ray interactions with the interstellar gas and from 
inverse Compton of ambient photons by cosmic ray electrons, and they have produced 
a map of the intensity distribution of the residual gamma-ray emission. Besides a few 
“point” sources, they hnd an excess in the central region extending somewhat North 
of the galactic plane, and a weaker emission from regions in the galactic halo. They 
mention an astrophysical interpretation for this halo emission: inverse Compton by 
cosmic ray electrons distributed on larger scales than those commonly discussed and 
with anomalously hard energy spectrum. 

I hnd it intriguing that the angular distribution of the halo emission resembles that 
expected from pair annihilation of dark matter WIMPs in the galactic halo (Gunn et 
al. 1978, Turner 1986), and moreover, that the gamma-ray intensity is similar to that 
expected from annihilations of a thermal relic with present mass density 0.1-0.2 of the 
critical density (Gondolo 1998a). Namely, the emission at 6 > 20° is approximately 
constant at a given angular distance from the galactic center - except for a region 
around (6, /) = (60°, 45°), correlated to the position of the Moon, and a region around 
{h,l) = (190°,—30°), where there is a local cloud {h and I are the galactic latitude 
and longitude, respectively). Dixon et al. (1998) argue against the possibility of 
WIMP annihilations on the base that direct annihilation of neutralinos into photons 
would give too low a gamma-ray signal. However most of the photons from WIMP 
annihilations are usually not produced directly but come from the decay of neutral 
pions generated in the particle cascades following annihilation. 
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Presently preferred dark matter candidates tend to give a gamma emission which 
is too low even in the continnnmP In a search for a snitable candidate it is worth 
examining the impact of various constraints on the properties of the candidate. I will 
later introduce a working model, so to make the discussion more concrete. 


2. A candidate WIMP can be a thermal relic. 


It is useful to use the WIMP mass and its annihilation cross section (times 
relative velocity at n = 0) av as parameters in the discussion. The requirement to 
approximately match the Dixon et ah maps selects a band in the av-m^ plane (band 
marked “halo 7 ’s” in £g. 1 ). Another band is selected by the requirement that the 
WIMP is a thermal relic from the early universe, with a relic abundance in the range 
0.025 < < 1 (band marked “0.025 < < 1” in fig. 1). The intersection of the 

two bands defines the interesting region. 

The gamma-ray band is obtained as follows. The gamma-ray intensity from WIMP 
annihilations in the galactic halo from direction with galactic longitude b and galactic 
latitude I is given by 

(j)^{bJ,>E) =n^{>E)^^ Jpldl. ( 1 ) 

X 

(l)^{b,l, >E) is in photons/(cm^ s sr), n^{>E) is the number of photons with energy 
above E generated per WIMP annihilation, av is the WIMP annihilation cross section 
times relative velocity, rriy. is the WIMP mass, and is the WIMP mass density in 
the halo. The integral in eq. (|I|) is along the line of sight in direction {b,l). This 
integral gives the angular dependence of the gamma-ray flux and depends on details 
of the dark halo model, which are not well known. For a canonical halo model, 
p(r) = pioc(?"c + where pioc is the WIMP mass density in the solar 

neighborhood, Vc is the halo core radius, R is the distance of the Sun from the galactic 
center, and r is the galactocentric distance. In this case, the integral depends only 
on the angle il) between the direction of observation and the galactic center. 
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Here x = 1 -I- {tc/R)'^ and c = cosTp = cos 5 cos Z. From dynamical studies one finds 
Pioc = 0.3-0. 5 GeV/cm^, R = 8-8.5 kpc, and Vc = 2-8 kpc. It is interesting to notice 
that the gamma-ray intensities aiip = 40° and '0 = 60° are approximately in the ratio 
2:1 as on the Dixon et ah maps. Assuming WIMP annihilation into quark-antiquark 
and lepton-antilepton pairs, fixing n.y(>lGeV) with the Lund Monte-Garlo, varying 

“Attempts to increase the flux by dumpiness in the halo (Berstrom, Edsjo and Ullio 1998) tend to 
produce either too many antiprotons or too few photons below 1 GeV, in contrast to the Dixon et 
al. maps (see section 4). 
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Fig. 1. Region in the annihilation cross section versus mass plane where the gamma-ray halo may 
be explained by annihilations of a thermal dark matter relic. 


the halo parameters in the range given above, and matching the observed intensity 
to eq. (0 to within 20%, I obtain the required WIMP annihilation cross section as a 
function of the WIMP mass. This is the band marked “halo y’s” in fig. 1. 

The second band in £g. 1 comes from the requirement that the WIMP relic density 
be in the cosmologically interesting range 0.025-1. {h is the Hubble constant in 
units of lOOkm/s/Mpc.) The WIMP relic density is related to the WIMP annihilation 
cross section through the approximate relation (Kolb & Turner 1990, Gondolo & 
Gelmini 1991) 


av 


2.0 X 10 ^^cm^/s 
ggV^XfVth^ 


( 3 ) 


where I have assumed that the annihilation cross section is dominated by the n = 0 
term both in the galactic halo and at freeze-out (s-wave annihilation). The freeze-out 
temperature Xfm can be obtained solving -b | ln( 5 f*/x/) = ‘&QA + \Ti{gm^av), with 
in GeV and av in cm^/s. For a Majorana WIMP g = 2, for a Dirac WIMP g = 4:. 
g^, is the effective number of relativistic degrees of freedom at freeze-out: ~ 81 
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before and ~ 16 after the QCD quark-hadron phase transition. Letting vary 
in the above range gives the relic density band in £g. 1. 

The two bands intersect for WIMP masses between 1.2 and 50 GeV. For example, 
G = 0.2 and H = 60 km/s/Mpc give the required cross section of 3 x 10“^® cm^/s at 
= 5 GeV. (This case was presented in Gondolo 1998a.) 


3. Constraints on candidates that couple to the Z boson. 


I assume in this section that yy annihilation and y-nucleon scattering are domi¬ 
nated by Z boson exchange. In this case, the annihilation cross section reads 
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where (3f = {1 — mj j, aj = Tf and Vf = Tf — 2efsm‘^ 9w are the usual axial 
and vector couplings of fermion / to the Z boson, and and are the analogous 
quantities for the y-Z coupling. 

The y-nucleon scattering cross section which is limited by negative direct 
dark matter searches, is related to the annihilation cross section by crossing sym¬ 
metry. The experimental bound on depends on the WIMP mass and on the 
spin-dependent or spin-independent character of the interaction (for a recent compi¬ 
lation of limits see Bernabei et al. 1998). For Z boson exchange, the spin-dependent 
and spin-independent y-nucleon cross sections read 
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Here = m^nij^/ {niy. + m^) is the reduced y-nucleon mass, and Aq is the quark 
q contribution to the spin of the proton. (From neutron and hyperon decay. An — 
Ad = 1.2573 ± 0.0028, while As is uncertain: As = 0 in the naive quark model. 
As = —0.11 ± 0.03 from deep inelastic data, and As = —0.15 ± 0.09 from elastic 
z/p —> z/p data.) 

The coupling of the y to the Z boson also gives a contribution to the Z boson 
decay width, if < mz- Namely, 


m ^ xx) 


CpTflz 

6V2n 


almlYl + Yx^l {ml + , 


(7) 


with Yx — (1 —4m^/m|)^/^. The experimental limit is V{Z xx) < 5 MeV (Barnett 
et al. 1996). 
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Once a relation between and is specified, the experimental bounds on 
and r(Z ^ xx) translate into a bound on av. Fig. 2 plots these bounds for two 
cases: a Dirac particle with {V — A interaction), and a Majorana particle 

with = 0 (axial interaction). For Dirac particles, the interesting region is not fully 
excluded by these constraints, but for Majorana particles it is. Hence the impact of 
these constraints is model dependent. 



Fig. 2. Constraints on the interesting region for a particle coupled to the Z boson. 


4. Constraints from dark matter searches. 

A candidate WIMP has to satisfy constraints from negative dark matter searches. 
In this section I consider indirect detection through production of rare cosmic rays 
(antiprotons) and through neutrino production in the Sun and the Earth, and direct 
detection through elastic scattering off nuclei in a laboratory detector. 

The most important constraint comes from the measured flux of cosmic ray an¬ 
tiprotons. If the gamma-rays are produced in jets originated by quarks, there is an 


5 





1 10 100 


[GeV] 

Fig. 3. Constraints on the interesting region from dark matter searches. 

associated production of antiprotons. The ratio of antiproton and gamma-ray fluxes 
is independent of the WIMP annihilation cross section and of the local mass den¬ 
sity, and the relative nnmber of antiprotons and photons per annihilation is hxed by 
the physics of jets. The antiproton flux at a p kinetic energy T at the top of the 
atmosphere is 

, diVp av 2 /o\ 

4>p{T) = Atiw? ( 8 ) 

X 

where dNp/dT is the antiproton spectrnm per annihilation, tcont is the p contain¬ 
ment time, which in the diffusion model of Chardonnet et ah (1996) is tcont — 
(1-|-p/3GeV)“°'®5 X 10^^ s, and fi = [T(T -|-2mp)]/[(T + A)(T -|-A-|-2mp)] takes into 
acconnt solar modulation. Fig. 3 shows the bonnd obtained by taking A = 600 MeV 
and imposing 0p(15O-3OOMeV) < 3 x 10“® p/cm^/s/sr/GeV (Moiseev et al. 1997). 
For WIMPs heavier than ~ 10 GeV, a related constraint on the gamma-ray flux be¬ 
tween 300 and 1000 MeV comes from the approximate relation n.y(300-1000MeV) ~ 
(iAp(150-300MeV)/(iT x lOGeV. This gives 0..y(3OO-lOOOMeV) < 1.3 x 10’^ pho- 
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tons/cm^/s/sr if photons originate in qnark jets with > 10 GeV. Compared with 
the observed value of ~ 8 x 10^ photons/cm^/s/sr, it implies too low a gamma-ray 
flux below 1 GeV. However, these bounds from antiprotons are uncertain because of 
uncertainties in the antiproton propagation in the galaxy and in the effect of the so¬ 
lar modulation. Moreover, these bounds depend sensitively on the relative branching 
ratios into the various annihilation channels, and for example can be avoided if the 
dominant decay channel is leptonic (see next section). 

The X fermions can also accumulate in the Sun and in the Earth, annihilate 
therein and produce GeV neutrinos. Accumulation is not efficient for WIMPs lighter 
than the evaporation mass, which is ~ 12 GeV for the Earth and ~ 3 GeV for the 
Sun. The curves marked “Earth” and “Sun” in Eg. 3 show approximate constraints 
obtained from the experimental upper bounds on the flux of through-going muons in 
the Baksan detector (Suvorova 1997). 

Another limit comes from direct searches through the crossing relation between 
annihilation and scattering cross sections. This relation can be written as 

(9) 

where c is the speed of light and / accounts for the details of the interactions and 
of the nuclear structure. The “direct detection” bound in Eg. 3 assumes / = 1, 
but actually / can range from 0 to inEnity. For example, consider a WIMP that 
couples to quarks and leptons only through exchange of scalar bosons, and consider 
annihilation through s-channel exchange and scattering through t-channel exchange. 
For kinematical reasons, if the scalar is GP-even, the annihilation cross section av at 
n = 0 vanishes and the scattering cross section at small velocity is Enite, so / = 0. On 
the other hand, if the scalar is GP-odd, it is the scattering cross section that vanishes 
and the annihilation cross section is Enite, so / is inEnite. Ghoosing the relative 
strength of GP-odd and GP-even exchange, the magnitude of the annihilation and 
scattering cross sections can be tuned. An explicit example is given in the next 
section. 


5. A candidate WIMP can be a sterile neutrino. 


The analysis in the previous paragraphs may have been instructive, but the con¬ 
clusions are so model-dependent that the discussion would be clearer if done in a 
speciEc model. Here I show that a particle with the required properties is a sterile 
right-handed neutrino Vg in a model with an extended Higgs sector. 

The Higgs sector contains two Higgs doublets Hi and H 2 and a Higgs singlet N. 
1 assume the following Higgs potential 


VEiggs = Ai [hIHi - vf) ' + A2 {hIH 2 - vt) ' + As 


N^-vl 
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A4 \H1H2 - ViV2\ + As H1H2 - V1V2 + -V 
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+ Ag 


hIh2 '. ( 10 ) 
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To fix the notation, H 1 H 2 = — H^H^- Taking all the A’s real and positive 

guarantees that the absolute minimum of the Higgs potential is at (H^) = Vi, (-H®) ~ 
V 2 , and (N) = vn, with the vacuum expectation values of all the other fields vanishing. 
Through the Yukawa terms 


■^Yukawa fdQ T fuQH2UR T feLH T , (H) 


Hi gives masses to the up-type quarks and the charged leptons, H 2 gives masses to 
the down-type quarks, and N gives a Majorana mass to the right-handed neutrino. 
There is no mixing of the right-handed neutrino with ordinary neutrinos, otherwise 
the new neutrino would have decayed in the early universe and would not be in the 
galactic halo at present.^ 

There are 2 would-be Goldstone bosons - a charged one and a neutral one 
G° - and 6 physical Higgs bosons: a charged one three neutral “scalars” Si 
(i = 1,..., 3), and two neutral “pseudoscalars” Pi {i = 1, 2). 

The physical charged Higgs field = H^* sin (3 + cos /? has squared mass 

m]j+ = Xqv"^. ( 12 ) 


As usual, V = \fvi + V2 and tan (3 = V2/V1. 

The squared mass matrix of the two physical pseudoscalar neutral Higgs helds 
reads 


/ (A4 -h As)^^ 2X5VVN \ 
y 2 X 5 VVN 4(A3 -h X4)v% J 


(13) 


in the basis v^(sin (3lmHi + cos ImiV). The orthogonal combination of fields 

is the would-be Goldstone boson G^. Let be the unitary matrix that diagonalizes 
namely Pi = U^A + U^Nj where i = 1,2. 

The scalar Higgs boson mass matrix is 


Ml = 


^ AXivl + (A4 -h A5)n2 
(A4 -|- X^)viV2 
V 2X5V2VN 


(A4 -|- X 5 )viV 2 
4 :X 2 V 2 + (A4 -h As)^! 
2X5VIVN 


2 X 5 V 2 VN 
2X5VIVN 
4(A3 -|- As)^^ 


(14) 


in the basis v^(ReiL°, ReH^, ReiV). Its mass eigenstates S'* (i = 1,..., 3) are obtained 
as Si = V2{UfiReHf + UgReH^ + UgReN). 

The original Higgs fields can be expressed in terms of the physical helds as 


Hi = Vi + i/l(ufiSi + iUfiPiSm/3), (15) 

H 2 = V 2 +^/l(U,^ 2 Si + iUfiPi cos f3) , (16) 

N = VN + ,/l{ugS, + tU,^Pi). (17) 

^F. Vissani has kindly pointed out that a particle that does not mix with ordinary neutrinos should 

not be called a “neutrino.” I keep this name for lack of a better one. 



This gives the interactions of the Higgs bosons with the quarks, the leptons and the 
right-handed neutrino, 


-^int — 


grriu, V 


2mw V]s[ ^ 
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_sin/3" 

grrid 

Ca < 

2mw 

COS f3 


Siuu + i cot iSU^Pi'u'y^u 
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It is then easy to work out the annihilation cross section. 


av = 


G^pirnl v"^ 
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TjP TjP 


TT vl ml^ - Ami 
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and the scattering cross section off nucleons. 




^uj\f 


TT (mjv + mu) v_ 
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\cos/5 siu/d^ 


n 2 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


Here c/ = 3 for quarks, c/ = 1 for leptons, Kf = cot (3 for up-type quarks, and Hf = 
tan f3 for down-type quarks and leptons. Moreover, kd = {mddd+msSS+mJ)h) = 0.21 
and ku = {mpau + m^cc + mitt) = 0.15. 

Table 1 lists important quantities for two interesting cases: a 3 GeV neutrino and 
a 7 GeV neutrino. 


Table 1. Two interesting cases. 



“3 GeV” 

“7 GeV” 

Experimental/ 

Observational 

Values 

mu [GeV] 

3 

7 


(/>7 (56°) [photons(>IGeV)/cm^ /s/sr] 

7.3 X 10-'^ 

7.4 X 10“^ 

~ 7.6 X 10-'^ 

H {H = 60 km/s/Mpc) 

0.11 

0.12 

~ 0.2 

(/>p(200MeV) [p/cm^/s/sr/GeV] 

~ 0 

7 X 10-® 

< 3 X 10-6 

O-XAA [Pb] 

1.2 X 10-3 

1.3 X 10“^ 

see text 

mp^,mp 2 [GeV] 

5,155 

4.2,156 


msi,rns 2 ,ms:i [GeV] 

5,155,220 

4.2,156,220 


mp± [GeV] 

no 

no 

>53 

T{Z^ffP,S)/T{Z^ff) 

1.3 X 10-® 

6.3 X 10“® 

<10-4 

“sin2(/3-a)” (Zh) 

2.3 X 10-^ 

2.3 X 10“® 

<10-2 

‘Wjd-a)” (hA) 

6 X 10-'^ 

6 X 10-® 

<0.3 


For the “1 GeV” model, the parameters are h = 0.9, v/vn = 52, tan/? = V 2 /V 1 = 
40, and Xi,..., \q = (the square of the SU(2) coupling constant). The gamma-ray 
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intensity is close to the observed one and the relic density is in the cosmologically 
interesting range. There is practically no antiproton flux from annihilations because 
only the channel is effective. The scattering cross section off nucleons is com¬ 

patible with the present bound of 2 x 10“^ pb. Constraints on the Higgs sector are 
discussed below. 

For the “7 GeV” model, the parameters are h = 0.8, v/vn = 20, tan/? = 33, 
and Ai,..., Ae = g^. The gamma-ray intensity is closed to the observed one, and the 
relic density is in the cosmologically interesting range. The antiproton flux is slightly 
higher than the measured one, but, given the big uncertainties in its estimation, it is 
compatible with it. The scattering cross section is smaller than the present limit of 
3 X 10“^ pb. The Higgs sector is discussed in the following. 

Two Higgs bosons {Si and Pi) are light and one has to worry about their pro¬ 
duction at colliders. But this production is very suppressed because the light Higgs 
bosons are mostly singlets. First, the LEP bounds on the search of Higgs bosons in 
two-doublet models are satisfied. Since there is an additional Higgs singlet, sin^(/5—a) 
in the bound from e+e“ — Zh must be replaced with |[/f^cos/5 -f P^siu/dp, and 
cos^(/5 — a) in the bound from hA must be replaced with — P/^Pfip. 

In the allowed region, the reinterpreted sin^(/5 — a) and cos^(/5 — a) are smaller than 
5x 10“^, and are not excluded by accelerator searches (Decamp et ah 1992). Secondly, 
Higgs bremsstrahlung from final state leptons and quarks in Z decays is suppressed 
by the small mixing between singlet and doublet Higgs bosons. In the present model. 


r(z ^ ffx) 

- //) 


x/2Gj. 



V™pA7. 


(23) 


where Ap^f = U^^cot/d, As^f = Pj^/cos/5 for up-type quarks, Ap.f = 

Asif = PjI/ siu/d for down-type quarks and leptons, c/ = 3 for quarks and c/ = 1 for 
leptons. The function g{y) comes from the phase-space integration, and is g{y) = 1 
at mx = 12 GeV. The values obtained in the two examples (see table 1) are lower 
than the LEP constraint on two leptons -|- two jets production r(Z f fX)/T{Z — 
//) < 1.5 X 10“^ at mx = 12 GeV (Decamp et ah 1992). 

I conclude that these two examples satisfy the experimental and observational 
constraints considered. A more general scan of the model parameter space is under 
way (Gondolo 1998b). 


6. Conclusions 

I have shown that it may be possible to explain the Dixon et al. gamma-ray 
emission from the galactic halo as due to halo WIMP annihilations. Not only the 
intensity and spatial pattern of the halo emission can be matched but also the relic 
density of the candidate WIMP can be in the cosmologically interesting domain. After 
a model-independent analysis to learn about the properties of a suitable candidate. 
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I have presented a working model: a sterile nentrino in a model with an extended 
Higgs sector. Two examples in parameter space indicate the existence of an interesting 
region in the model parameter space in which present observational and experimental 
constraints are satisfied and the gamma-ray emission is reprodnced. 
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